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ABSTRACT 

We study the formation of molecular hydrogen in cooling gas behind shocks 
produced during the blow-away process thought to occur in the first collapsed, 
luminous (Pop III) objects in the early universe. We find that for a wide range 
of physical parameters the H 2 fraction is / 6 x 10~ 3 . The H 2 mass produced 
in such explosions can exceed the amount of relic H 2 destroyed inside the 
photodissociation region surrounding a given Pop III. We conclude that these 
first objects, differently from the suggestion of Haiman et al. 1997, might have a 
net positive feedback on subsequent galactic formation. We discuss the effects of 
radiation and the implications of our results for the soft-UV background. 

Subject headings: Cosmology: theory - galaxy formation - molecules - 
supernovae 
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1. Introduction 

Current models of cosmic structure formation based on CDM scenarios predict that the 
first collapsed, luminous (hereafter Pop III) objects should form at redshift z 30 and have 
a total mass M « 1O 6 M or baryonic mass M b & 1O 5 M (Couchman & Rees 1986, Haiman 
et al. 1997 [HRL], Tegmark et al. 1997). This conclusion is reached by requiring that the 
cooling time, t c , of the gas is shorter than the Hubble time, tn, at the formation epoch. 
In a plasma of primordial composition the only efficient coolant in the temperature range 
T < 10 4 K, the typical virial temperature of Pop III dark matter halos, is represented by 
H 2 molecules whose abundance increases from its initial post-recombination relic value to 
higher values during the halo collapse phase. It is therefore crucial to determine the cosmic 
evolution of such species in the early universe to clarify if small structures can continue to 
collapse according to the postulated hierarchical structure growth or if, lacking a cooling 
source, the mass build-up sequence comes to a temporary halt. 

The appearance of Pop III objects is now thought to cause a partial destruction of the 
available molecular hydrogen either in the intergalactic medium (IGM) and/or in collapsing 
structures; the result is a negative feedback on galaxy formation. This effect has been 
pointed out by HRL and it works as follows. As stars form in the very first generation of 
objects, the emitted photons in the energy band 11.2-13.6 eV are able to penetrate the gas 
and photodissociate H 2 molecules both in the IGM and in the nearest collapsing structures, 
if they can propagate that far from their source. This negative feedback and its possible 
limitations are discussed by Ciardi, Abel & Ferrara (1998, CAF). Here we propose and 
investigate instead a possible positive feedback based on supernova (SN) explosions, which, 
under many aspects, is reminiscent of a scaled version of the explosive galaxy formation 
scenario introduced by Ostriker & Cowie (1981) and put forward by many others. Pop Ills 
are very fragile due to their low mass and shallow gravitational potential: only a few SNe 
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are sufficient to blow- away (Ferrara 1998) their baryonic content and drive an expanding 
blastwave into the IGM, which eventually becomes radiative and allows the swept gas to 
cool in a dense shell. The cooling transient, as we will see, is characterized by a strong 
nonequilibrium condition in which recombination lags behind the temperature decrease. 
As already pointed out by Shapiro & Kang (1987) and Kang & Shapiro (1992), this is a 
favorable condition for H 2 formation. Our conclusions are that the amount of molecular 
hydrogen thus formed can exceed the one destroyed via photodissociation, yielding a net 
increase of the H 2 in the universe. As a consequence, not only the galaxy formation process 
is not halted, but instead it is favored by the effects of Pop III formation. 

Sec. 2 describes the main properties of multi-SN explosions in the early universe; Sec. 
3 is devoted to the calculation of H 2 formation in their cooling shell. In Sec. 4 we compare 
the magnitude of positive and negative feedbacks; some discussion of the results is given in 
Sec. 5. 



2. Properties of Pop III Remnants 

The mechanical luminosity of an OB association in a Pop III object at redshift z can 
be written as (Ciardi & Ferrara 1997, CF) 

L = ~ 4 x 10 3f U, 5 / b , 8 (l + ^ 2 M 6 ergs s" 1 ; (1) 

Tt ff 

where (1 + -2)30 = (l + z)/30 and e = 10 51 ergs is the energy of a SN explosion; we assume a 
Salpeter IMF, according to which one supernova is produced for each 56 M = z/ -1 of stars 
formed. The baryon density parameter is Q b = 0.05f2{, i5 , of which a fraction f b ~ 0.08/^g 
(Abel et al. 1997b) is able to cool and become available to form stars. The dark matter halo 
mass and density are M 6 = M/1O 6 M and p ~ 200p c = 200[1.88 x l(T 29 /i 2 (l + zf] g cm" 3 , 
respectively (we assume a flat EdS universe Q m = 1 and h = H/100 km s _1 Mpc -1 ); 
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the corresponding free-fall time is tff = (AnGp) -1 ! 2 ; r _1 = 0.6% is the star formation 
efficiency, calibrated on the Milky Way The correlated SN explosions drive a blastwave 
in the surrounding gas, which eventually propagates into the IGM; due to their low mass, 
the effects of Pop III ISM on the blastwave expansion can be neglected (MacLow & 
Ferrara 1998). The evolution of the shock radius, R s , can be obtained in the thin shell 
approximation by solving numerically eqs. (5)- (6) of CF. For sake of simplicity we use the 
following (Sedov) analytical solution of the same equations which holds for the adiabatic 
case when the external confining pressure and the gravitational pull of the DM halo are 
neglected: 

{Et 2 \ 1/5 

R s (t) ~ 0.76 - , (2) 

where E = Lt OB , t OB ~ 10 7 yr is the average lifetime of massive stars, E is the total energy 
of the explosion, p^ = Q^p the IGM density at redshift z. Eq. |2| does not take into account 
the effects of Hubble expansion, which are negligible since R s <C c/H (see eq. |]). 

The corresponding temperature of the postshock gas is 27.6(i? s /kms -1 ) 2 K. Such hot 
gas will not be able to cool until its cooling time, i c , becomes smaller than the Hubble 
time, iff. The energy losses, A(T), due to gas-radiative processes in a primordial gas are of 
order 10~ 23 ergs cm 3 s -1 in the range 10 5 < T < 10 8 K (Schmutzler & Tscharnuter 1993); 
the other relevant cooling agent is inverse Compton cooling off CMB photons (Ikeuchi & 
Ostriker 1986). Comparing the two rates, it is easily seen that inverse Compton dominates 
the cooling for z ^ \85Qbh 2 — 1, for T = 10 5 K. The condition i c < iff requires that 
z £ 8/i 2 / 5 — 1; shocked gas produced at epochs earlier than this will be able to cool and 
condense in a dense shell behind the shock front. Note that as the temperature drops to 
lower values with subsequent gas recombination, cooling via hydrogen Lya line excitation 
will start to become important. At i = i c = 1.2 x 10 3 (1 + z)~ 4 Gyr, the shock radius and 
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postshock gas temperature are respectively 

R s (t c ) = 0.16/^ 5 (l + ^ 3 - 19/10 M 6 1/5 ^ 1/5 kpc; (3) 

T s (t c ) = 4.9 x 10 4 /S 5 (1 + z)H /5 M 2 6 /5 h- 2 ^K. (4) 

These relations hold if the Mach number of the shock with respect to the IGM sound speed 
is M. 3> 1, a condition verified if the IGM has not yet been reheated. 

3. Molecular Hydrogen Formation 

In the following, we calculate the time-dependent, nonequilibrium ionization and 
molecule abundance evolution behind the multi-SN shock as a function of the object 
mass and epoch of the explosion. The postshock cooling gas is assumed to be initially at 
temperature T s (t c ) and density 4p fe . In hydrogen of density n = nn + n H + + n# 2 (we neglect 
He), the evolution of the ionization fraction, x = nn+/n, molecular fraction, / = nn 2 /n, 
and temperature, T, is determined by the following rate equations: 

x = kinx(l — x) — k 2 nx 2 (5) 



/ = k 8 n H -(l - x) - k u nxf - [k 12 x + k 13 (l - x)]nf (6) 



( 7 -l)A(T) x 
T = -T ^ i^^. + . 7 

7 p 1 + x 

The various rate coefficients ki are labelled according to the nomenclature given in Abel et 
al. 1997a (Appendix A); n H - is the H~ species number density; 7 = 5/3 is the specific heat 
ratio and p — (1 + x)nksT is the gas pressure. In addition to the standard atomic cooling 
processes, A(T) includes Compton and H 2 cooling; the latter has been adopted from Martin 
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et al. (1996). We assume that no heating is provided to the gas after it has been heated by 
the shock to the initial temperature and neglect adiabatic cooling due to Hubble expansion 
for reasons similar to those given in the previous Section. 

The first equation describes the usual collisional ionization (rate k\) - recombination 
(k 2 ) balance. The equation for / is slightly more complicated and we describe it shortly. 
At the low densities considered here, in the absence of dust, H 2 is formed in the gas phase 
mainly via the channel H + e~ — > H~ + hv, at rate h$ (formation via the H 2 channel, when 
included, is found to be negligible in our case); it is then destroyed by one of the following 
mechanisms: (i) it is collisionally ionized by proton impacts (k u ); (ii) it is dissociated 
by collisions with electrons (k 12 ) and/or H-atoms (k 13 ). Abel et al. 1997a have shown 
that using the equilibrium value for nn- is a very good approximation since the reactions 
determining the H~ abundance occur on a much shorter timescale than those responsible 
for the if 2 chemistry; therefore we have calculated n H - from their eq. (24). Also, we 
note that in the above form, the rate equations correspond to the so-called minimal model 
described in that paper. The energy equation has been solved in the isobaric limit; this 
approximation is known to be appropriate for the evolution of Galactic supernova remnants 
(Cioffi et al. 1988), where most of the radiative phase is spent at constant pressure; the 
same result was found to hold by Shapiro & Kang (1987) and Kang & Shapiro (1992) for 
a wide range of shock velocities in a cosmo logical context. The numerical solution of the 
above rate equations are shown in Fig. 1 for two different object masses (M 6 = 1, 10) 
and initial redshift (z = 20,25,30). Almost independently on the specific values of these 
parameters a high H 2 fraction (/ ~ 6 x 10~ 3 ) is produced by redshift ten, with a steep 
increase occurring when the gas has cooled down at T 20, 000 K. H 2 production is 
faster in high-z, high T s shocks, but the peak abundance is practically the same in all 
cases; the final gas temperature is in the range 300 — 500 K. We have also experimented 
a different H 2 cooling function (Hollenbach & McKee 1979): this produces a cooler final 



- 8- 



state (T = 80 — 200 K) but only a marginally different asymptotic value /si5x 10~ 3 . For 
comparison, we have also plotted in Fig. 1 the Tegmark et al. 1997 "rule of thumb" value 
/ = 5 x 10~ 4 at which H 2 cooling becomes important. 

4. Positive or Negative Feedback ? 

The results obtained in Sec. 3 show that typically at least a H 2 fraction / 6 = //6 x 10 -3 
is formed after explosive events leading to the blow-away of Pop Ills. As discussed in the 
Introduction, HRL have argued that the same objects could suppress the surrounding H 2 
abundance due to their UV photodissociating radiation. In order to evaluate the impact 
of Pop Ills on the subsequent galaxy formation, largely regulated by the availability of 
H 2 in this mass (and redshift) range, it is useful to compare the H 2 mass production vs. 
destruction. A lower limit to the amount of H 2 produced in an explosion is readily found 
to be equal to M^(z) ~ (ATr/3)p b R s (t c ) 3 (2f) or 

M+ (z) ~ 102/g 5 fi b , 5 (l + z)3 O 27/lo M 6 1/5 /6^ /1O M . (8) 

As shown by CAF, the UV/ionizing radiation from Pop Ill's massive stars previous to 
blow-away will produce both a HII region and a region of photodissociated intergalactic H 2 
(radius Rd) in which the object is embedded. The radius Rd can be defined by requiring 
that the photodissociation timescale (td ~ k 2 j , notation as in eqs. ||||) is shorter than tjj- 
This condition yields the definition R d = S]^(l + z)~ z ^h~ 1 / 2 , where Slw — /3Si(0) is the 
UV photon flux in the H 2 Lyman-Werner (LW) bands (11.2-13.6 eV), assumed here to be 
proportional to the flux of Lyc photons, Si(0), just before the massive stars explode. The 
value of the constant (3 depends somewhat on the IMF and on the evolutionary stage of 
the stellar cluster, but its value should be reasonably close to unity. Paralleling eq. [I], we 
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estimate S^O) to be 

m = fu VpP fesM b h M ^ S-\ (9) 

Tt OB m p 

where f uvpp = 48000/ m , PP) 48 is the UV photon production per collapsed proton efficiency 
(Tegmark et al. 1994) and f esc = 0.2f esC:20 is the escape fraction of such photons (absorption 
can be caused both by interstellar neutral H and/or dust). We have checked that this simple 
estimate is within a factor of 2 of the value obtained from the recently revised version of the 
Bruzual & Chariot (1993) spectrophotometric code. The adopted reference value f esc = 0.2 
is an upper limit derived from observational (Leitherer et al. 1995, Hurwitz et al. 1997) and 
theoretical (Dove & Shull 1994) studies. It follows that 

R d = 2.4( / 3S 48 ) 1 / 2 (l + z)-^' 1 / 2 kpc, (10) 

where S^s = 5'j(0)/10 48 s _1 . Thus, the ratio between the H 2 mass produced and destroyed 
by Pop Ills is 

MU Z ) = ( f ) 
Mn 2 (z) \ fiGM ) 

The post-recombination relic fraction of intergalactic H 2 is estimated to be 
fiGM ~ 2 x 10~ 6 /i -1 (Palla et al. 1995; Anninos & Norman 1996; Lepp, private 
communication). The previous relation is graphically displayed in Fig. 2, along with the 
values of R s (t c ) and Rd- From that plot we see that objects of total mass M 6 = 1 produce 
more H 2 than they destroy for z £ 25; larger objects (M 6 = 10) provide a similar positive 
feedback only for z £j 15, since they are characterized by a higher Rd/R s (t c ) ratio. However, 
since in a hierarchical model larger masses form later, even for these objects the overall 
effect should be a net H 2 production. 



R s (t c 
Rd 



(11) 
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5. Discussion and Conclusions 

The results obtained in the previous Section could be modified by the fact that so far 
we have neglected the effects of the Pop III stellar cluster radiation impinging on the cooled 
shell. In fact, even after all the SNe have exploded blowing-away the gas, the coeval low 
mass stars will continue to produce some residual flux in the energy range 0.755 — 13.6 eV, 
relevant to the H 2 formation network. Two processes might have some effect in this context: 
(i) H~ photoionization, and (ii) H 2 two-step photodissociation. 

The first process occurs at a rate k 23 ~ 6 x 10~ 17 Siw/R 2 (t c ) s" 1 per H~ atom, or 
2.4 x 10- 10 (3Mq /5 s" 1 using eqs. [3], |9| and reference values for the remaining parameters. 
The most effective H~ destruction process among the ones included is H 2 formation, whose 
rate in the same case is ~ 1.4 x 10~ 9 s -1 . Thus, photoionization is negligible, and we have 
checked this conclusion numerically. 

We can calculate the effects of two-step photodissociation on the shell as follows. 
The shell H 2 column density is N Ha ^ f M s /m p R s {t c ) 2 = 1.5 x 10 17 (1 + z)H /W Mq /5 cm" 2 
(we have assumed that all the swept IGM mass is in the cool shell). Since the critical 
column density for self-shielding is N^ 2 lt = 10 14 cm -2 (Draine & Bertoldi 1997), the 
shell is optically thick to dissociating radiation. Nevertheless, the flux in the LW bands, 
Jd = SLwh p /4:iTRg, where h p is the Planck constant, will induce a dissociation front in the 
shell; Jd is very likely dominating on a possible background radiation. The propagation 
speed will be ~ £\d/td, where = /n H2 is the LW photon mean free path, and 
£ ~ 3 — 5 is a constant obtained by comparison with numerical results. Thus, the time 
required to dissociate the entire shell is ~ NH 2 /nH 2 Vd = 3 x 10 8 (Jd,2iO~ 1 -^H 2 ,i7 J 1 ^ where 
<-7d,2i = </d/10~ 21 ergs cm 2 s~ 1 Hz~ 1 sr" 1 , and Nh 2 ,x7 = N H2 /10 17 cm~ 2 . Since theoretical work 
(Bruzual & Chariot 1993) suggests that /3 drops to ~ 3 x 10~ 3 <C 1 after all the massive 
stars (the main contributors to the UV flux) have died, dissociating the entire shell will 
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require a very long time, certainly larger than the Hubble time. Finally, Kang & Shapiro 
(1992) pointed out that any external ionizing radiation field (which could be provided by 
the hot postshock gas and/or by the residual intermediate mass stars in the cluster) tends 
to increase the final value of /, although introducing a time delay due to a temperature 
plateau in the evolution characteristic of photoionization heating. Thus, our estimate seems 
to provide a robust lower limit to the molecular hydrogen abundance. 

How can this H 2 enhanced gas be used for galaxy formation ? First, regions in which 
Pop Ills are clustered enough for their shells to interact might become sites of active 
star/galaxy formation. Shell interactions appear to be necessary to make them gravitational 
unstable, since the temperature required for isolated fragmentation (Couchman & Rees 
1986) is well below the values ( £ 100 K) found here. Next, these intergalactic shells can be 
accreted by neighbor objects, possibly increasing their H 2 abundance above the threshold 
/^5x 10 -4 required for efficient cooling (Tegmark et al. 1997). 

We conclude that Pop III objects can produce regions of considerably high molecular 
hydrogen abundance due to multi-SN shocks propagating in the IGM. We have also 
seen that the H 2 thus produced can exceed the amount of relic H 2 destroyed inside 
the photodissociation region surrounding a given Pop III. This occurrence suggests that 
these first objects might have a positive feedback on subsequent galactic formation. Such 
conclusion is certainly valid on a local scale, as defined by the radii of the two influence 
spheres, R s and Rd, of a single object. However, our definition of Rd (see Sec. 4) does not 
exclude that some photodissociating flux is present outside Rd- If many of these objects 
are present in the universe at the same epoch, they might contribute to an early soft-UV 
background which, illuminating an isolated, collapsing cloud, could photodissociate and 
depress its H 2 content; this possibility has been discussed by HRL. However, CAF have 
demonstrated that, once the surviving relic IGM fraction between photodissociated spheres is 
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taken into account, the universe opacity to LW photons before reionization becomes of order 
unity on scales larger than R d , but smaller than the typical interdistance between Pop Ills 
at z = 20 — 30. This implies that the soft-UV background is very weak due to intergalactic 
absorption, and, as a consequence, the proposed negative feedback considerably inhibited. 
Furthermore, most of the dissociating flux is absorbed by the thick H 2 shells discussed 
above. Given that in general Rd > Ri, where Ri is the HII region radius of a Pop III 
(CAF), clearing of intergalactic H 2 by photodissociated sphere overlapping will occur before 
complete reionization. It is not clear at this stage if that event occurs before the typical 
mass scale of the collapsing objects is such that cooling can proceed via Lya only, thus 
considerably weakening the arguments in favor of a temporary halt of galaxy formation at 
high redshift. 



We acknowledge T. Abel, B. Ciardi, D. Cox and G. Field for discussions; E. Corbelli, 
Z. Haiman, A. Loeb, F. Palla, O. Pantano, P. Pietrini, and Y. Shchekinov for comments; 
thanks are due to L. Pozzetti for help with the spectrophotometric code. 
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redshift 



Fig. 1. — Evolution of the H 2 molecular fraction / (solid), H ionization fraction x (dotted) 
and temperature T/10 6 K (dashed) in the postshock gas for different values of the explosion 
redshift z = 20, 25, 30 and Pop III mass M = 1O 6 M (upper panel) and 1O 7 M (bottom). 
The cosmological parameters are fi& = 0.05, h — 1; the cooled baryon fraction is fb = 0.08. 
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Fig. 2. — Ratio between the H 2 mass formed and destroyed by a Pop III object as a function 
of the multi-SN explosion redshift (solid lines); values larger than one for the ratio define 
the epochs where Pop Ills have a positive feedback on galaxy formation. Also shown are the 
shell (proper) radius at cooling, R s (t c ) (dotted), and the photodissociation (proper) radius, 



